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The behavior of a Josephson flux-flow oscillator in the presence of both bias current 
and magnetic field fluctuations has been studied. To derive the equation for slow phase 
dynamics in the limit of small noise intensity the Poincare method has been used. Both 
the form of spectral line and the linewidth of the flux-flow oscillator have been derived 
exactly on the basis of technique presented in the book of Malakhov |^] jl|, known limiting 
cases are considered, limits of their applicability are discussed and appearance of excess 
noise is explained. Good coincidence of theoretical description with experimental results 
has been demonstrated. 
PACS number(s): 74.50. + r, 7 A AO. + k 

I. INTRODUCTION 

Long Josephson oscillators operating in the flux-flow regime O are presently considered as possible 
devices for applications in superconducting millimeter- wave electronics In comparison to single fluxon 
oscillators they have higher output power, wider bandwidth, and easier tunability, but they have a wider 
linewidth Q of the emitted radiation from the junction. Recent measurements by Koshelets et al. |||| 
have indeed shown a linewidth for a Josephson flux-flow oscillator which is of about one order of magnitude 
wider than the one derived for a short (lumped) Josephson junction This last property is quite 

undesirable if one wishes to use such devices, for example, as local oscillators in radioastronomy receivers 
||. For concrete applications it is important to get a model which adequately describes the linewidth 
of the flux- flow oscillator. With such a model one can hope to control the phenomenon of linewidth 
broadening by properly choosing the design parameters of the device. A first attempt in this direction 
was performed by Golubov et al. and Ustinov et al. [|l0|Jll| 1 in terms of a particle model for the train 
of fluxons moving in the junction. However, the authors of JTof calculated the variance of frequency 
fluctuations (that is quantity difficult to measure) but not a linewidth. Another attempt to derive the 
linewidth was performed in |l2| , but the results obtained are restricted by the consideration of " particle- 
like" picture of fluxon motion in an infinite junction and the linewidth is expressed in quantities that 
are not easily accessible from experiment in a flux-flow regime, e.g., average interval between fluxons. 
Besides, magnetic field fluctuations and parametric effects that may lead to additional broadening of 
the linewidth, were not considered in those papers. Recently, importance of accounting magnetic field 
fluctuations has been discussed in [ fi"3|| . 

The task of deriving of the linewidth of FFO may be decomposed into two parts: one is more experi- 
mental and another one is more methodological. 

One of the difficulties of the considered problem is absence of understanding of nature of noises of FFO. 
It is clear, that there are several different noise sources influencing the FFO: natural wideband noises 
(such as thermal and shot noises), technical narrowband noises and possibly flicker noise. And all these 
noise sources affect the FFO via fluctuations of both bias current and magnetic field (control line current). 
Considering present FFO designs |l4f| , one can guess, that there are some noise components in bias and 
control line currents that are correlated and some that are uncorrelated and lots of detailed experimental 
study should be performed to understand nature of this complicated mixture of fluctuations. 

On the other hand, if we know the parameters of noise (say, parameters of both natural and technical 
fluctuations), our task is to obtain the linewidth and the form of spectral line of FFO, and on the basis 
of the obtained characteristics to predict how to improve the noise properties of the oscillator. And here 
we would like to consider this "methodological" part of the analysis of noise properties of FFO. 

The aim of the present paper is to give strictly mathematical derivation of the required fluctuational 
equation for slow component of the phase that can be done for the most important case of small noise 
intensity, and present exact derivation of the linewidth and the form of spectral line on the basis of 



1 This paper is dedicated to memory of my teacher, The Honored Scientist of Russia, Russian State Prize Winner, 
Prof. Askold N. Malakhov (5.12.1926-7.11.2000). 
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methods described in the book by Malakhov We will assume a certain model of noise sources of 
FFO bias and control line current fluctuations, as well as we will consider the parametric effect of higher 
harmonics leading to additional broadening of the spectral line. 



II. BASIC EQUATIONS 

The electrodynamics of a long Josephson junction in the presence of magnetic field is described by the 
perturbed sine-Gordon equation 



subject to the boundary conditions 

dc/)(0,t) _ dcj)(L,t) 
dx dx 



= r. (2) 



In this equation space and time have been normalized to the Josephson penetration length Xj and to the 
inverse plasma frequency w" 1 , respectively, a is the loss parameter, r\ is the normalized dc bias current 
density and T is the normalized magnetic field. In accordance with RSJ model |j| one takes the loss 
parameter a = -f-, where uj p = y^2eI c /hC, lu c — 2eI c Rpf /h, C is the capacitance, Rpj is the normal state 

resistance (Rn — V/I qp , V being voltage and I qp - the quasiparticle component of the current), I c is the 

i i 

critical current, r\ = J / J c {I = J J{x)dx, I c = J J c (x)dx, I is the bias current), I is dimensional length 

o o 

of the junction, L = l/Xj. 

In general, both bias current r\ and magnetic field T (control line current) may fluctuate: rj = rio + 
r]F{x,t), r = To + Tp(x,t) and usually these fluctuations are supposed to be wideband noises and are 
small: for 770 ^ and Tq ^ variances of rjp{x, t)/r]o, Tp(x, t) /Tq are much smaller than unity. Therefore, 
we will consider the noise sources rjF(x,t) and Tp(x,t) as perturbations that do not affect the current- 
voltage characteristic, but lead to nonzero width of the spectral line. Following ||, we suppose that 
r]F{x,t) is Gaussian noise with zero mean value (r]F(x,t)) = and its spectral density is so wide that 
r]F{x,t) may be treated as white noise with the correlation function: 

( VF (x, t) V F(x', t')) = ™f"l 8(x - x')S(t - t'). (3) 

Here and in the following <> denotes ensemble average, ks is the Boltzmann constant and T is the 
temperature. In comparison with j7|, we consider simple RSJ model for current fluctuations: usually, at 
standard working temperature T = A.2K , pair current fluctuations are much smaller than quasiparticle- 
current fluctuations and may be neglected: I p = 0. Also we do not take into consideration the shot noise 
contribution that may be neglected if the condition 2ksT ^ eV is fulfilled. 

The properties of the magnetic field fluctuations TF(x,t) were not studied in the literature. From 
the present designs of FFO §,|,@ one can, however, make some conclusions about nature of these 
fluctuations. In the present layouts the base electrode of the long Josephson junction is employed as a 
control line. Therefore, wideband bias current fluctuations will enter the control line. Moreover, following 
recent idea of Koshelets (experimentally confirmed in [Q), even if the control line is isolated from the 
junction, fluctuating bias current may induce magnetic field, that will affect fluxons. On the other hand, 
narrowband technical fluctuations also exist there. So, we can model the control line fluctuations as 
follows: Tp(x,t) = arip{x,t) + Ti(x,t) + Tx{x,t), where Ti(x,t) are internal control line pair-current 
fluctuations (we neglect by Tj(x,t) = 0, supposing that they are much smaller than ar]F{x,t)) and 
Tt{x, t) are narrowband technical fluctuations. Since there are many compensation techniques that allow 
to significantly eliminate influence of narrowband technical fluctuations |],[|J|[l3|| , we will also neglect 
them: Tx{x,t) — 0. The question about attenuation factor a via which bias current fluctuations are 
converted into magnetic field fluctuations is not trivial and lot of theoretical and especially experimental 
work should be done to answer this question. It is clear, that a will be different for different types of 
FFOs and depends on the junction geometry and distribution of currents in the base electrode. Let us 
suppose, that the value of a is known and later we will discuss how a may be measured. We note, that 
the approach for linewidth calculation |l[] presented below recently was successfully used for calculation 
of the linewidth of Cherenkov FFO pi . This approach is rather universal and allows to take into account 
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almost any noise sources, even nicker noise (for which power spectral density diverges for lu — > 0, but, 
nevertheless, calculation of the linewidth may be done in this case [jjj) and in the present paper we neglect 
technical fluctuations only to simplify the analysis. 

The flux-flow regime is characterized by excitations which travel on top of a fast rotating background so 
that the effective nonlinearity in the system is drastically reduced due to fulfilling the following conditions: 
rj/a 3> 1 and (or) r > 1. In order to derive the linearized equation for slow component of the phase 
<p(x, t) (that is required for obtaining the spectral characteristics) we will first derive it in the case of zero 
noise intensity (j]F(x,t) = Tp(x,t) = 0) and later will consider noise as small perturbation. In papers 
Pq| , [Til linear mode theory and perturbative analysis around rotating background ((f) = (j>o + tp, ip <C 1) 
have been used to derive the current- voltage characteristic of FFO. 

We will use more general Poincarc method: obtain the solution as the series with respect to the small 

parameter e = ( f )' « 1. Let us change variables in Ec, §), r=^ Z = fr. 

d 2 (/> d(f> d 2 (j) . 

d^ +(3 ^-d^^ P - £Smi ^ (4) 

where (3 — a 2 /r]. 

The steady-state solution of this equation we will find in the form: (/>(r) = 4>o(t) + e<fii (r) + £ 2 </>2(t) + . . . 
(|0o (t) I 3> e l0i( T )| ^ e 2 1 02 ("7") | ^ . . .). Substituting this into Eq. (||) we will find the zero order equation: 

<9 2 o d(j) Q d 2 (j) 

~d^ + ^-^= P - (5) 

It is easy to see, that the steady-state solution of this equation is: 4>o(t) — t + jz — + Tx, 7 = aT/rj. 
To get higher order equations we have to decompose sin(0o(T, z) + e0i(r, z) + e 2 02(T, z) + . . .) into Taylor 
expansion. From the structure of the considered linear recurrent equations we know, that the steady-state 
solution 0„(r, z) may be presented in the form: n (r, z) — uj n r + (f> np (T,z), where </> np (r, z) is periodic 
nongrowing component. 

Let us now collect together all linearly growing components co n T and we will get: sm(0(r, z)) = 
sin({woT + euj\T + e 2 uj2T + . . . + 72} + e4>i p (T) + e 2 (f>2p{'r) + ...). Now we can linearize sin(</>) as: sin(0) w 
sin(w J T+-fz) + e(f) 1 p(T, z) cos(uijt + jz) + e 2 2 p(r, z) cos(ujjT + -fz) + . . ., where ujj = ui + euji + e 2 uj2 + ■ ■ ■ is 
the oscillation frequency (loq — I), and 10%,..., oj n ,..., <Pi p (t, z), 02 P (t, z), </>„ p (t, z), ... are unknown 
functions that we want to obtain. Restricting ourselves by consideration the solution up to the 2-nd 
order only (in principle we can do it up to any order, all equations may be solved recursively), we get the 
following equations for 0i(t, z), fair, z): 

d 2 <j>i , R d(t>i <9 2 </>i ■ / . \ 

-g^2-+^--^2- = -^ip(r)cos(cjjT + 1 z). (7) 
It is easy to see from Eq. (^|) that u>i — and substituting the solution in the form 

oo 

4>i p (t, z) = [Ai n cos(cjjt) + Bi n sin(wjT)] cos(k n z) 

n=0 

into (|^) one can find Ai n and Bi n : 

Xi n = {2-6 , n ) ^ jIc -&- Jij)Is , (8) 

Bi n = -^-S J UJjIs + ^- U}2j)Ic , (9) 

{^jy + {k n -u 2 y 

1 C 1 f^ 1 

Is = = sin (72;) cos(k n z)dz, Ic = = / cos(7z) coa(k n z)dz, L = — L. (10) 
L Jo L Jo a 
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Substituting this solution (where uj = 1 + e 2 uj 2 ) into Eq. (Q), we get: 



w 2 



E 



2 — #0,n 



(l + e 2 U, 2 )[J§+I§] 



.—2 



(0(1 + e 2 ^)) 2 + [fc„ - (1 + e 2 W2 ) 



(11) 



From this equation k>2 may be found. Analogically to (f>i p (t) one can find <j) 2p (t) . 

Combining Eqs. (§|)-(@) together we can get equation for ip(r, z) — <po(r, z) + e</>i(r, z) + e 2 <fi 2 (T, z) + 
. . . + e"0 n (r, z) that is equivalent to the steady-state case of the original equation (Q) up to the n-th order 
and lim ip(r, z) — <P(t, z): 



d 2 ib dib d 2 ib 
'■ft 



(3 — e sin(wjT + jz) 



2L 



p CO 

/ pin cos(7z) - Bi n sin(7z)] cos(k n z)dz 

n n=0 



2 00 



^ pin cos(2wjr + 72) + Bi n sm(2oj jt + 7z)] cos(knz) - . . . , (12) 



n=0 



where ojj = 1 + ewi + e 2 o;2 + . . . + e n ui n is supposed to be known. Here we have presented in the explicit 

form the term 4>i p (t) cos(wjt-\-^z) and have taken into account that only zero mode term with m = will 

00 

contribute into linearly growing component of 4>%(t) = ^ A 2 „ l (t) cos(k m z) , whereas all contributions 



m— 



L 00 



due to = / Pin cos(72;) — Bi n sin(7z)] cos(fc„z) cos(k m z)dz, m ^ 0, will decay with time. 



n=0 



Now, before introducing noise sources, it is convenient to change variables back: 

d 2 ifj dip d 2 ip 



dt 2 



where f(x,t) = f s {x,t) + f c (x,t), f s (x,t) = sin(Cljt + Tx), r\ + afl 2 = aClj, 

^ :v 

f c (x, t) = - [Ain cos(2fiji + Tx) + Bin sin(2Qjt + Tx)] cos(fc„a 



n=0 



(13) 



(14) 



1 f 1 00 

afi 2 = — / - [^4in cos(rx) — £?i n sin(ra;)] cos(k n x)da 



(15) 



^4i„ and Bin are: 



^4ln = (2 - S , n ) 



aUjIg - (fc 2 - n 2 j)I s 

( a nj) 2 + (k 2 n-n 2 ) 2 ' 



(16) 



o __/ 9 _r , anjJs + (fc^-^)/ c 

ln ~ 1 (aa 7 ) 2 + (fc 2 - tt 2 ) 2 ' 



FL sin(TL) cos(7rn) 



(ri) 2 -(7rn) 



2 — , *S 



rL(l-cos(rL) cos(Trn)) 



(17) 



(rL) 2 -(7rn) 2 



fcn = im/L, CI j — rjcjj/a is the oscillation fre- 



and Ic 
quency. 

Let us focus on the equation for the second-order frequency correction (|l l|) , that in original variables, 
substituting explicit form of Ic and Is, looks like: 



(2 - 6 . n )(ri/a + fl 2 )(TL) 2 [l - cos(TL) cos(Trn)] 
^ [(r, + an 2 ) 2 + (fc 2 - ( V /a + n 2 ) 2 ) 2 }[(TL) 2 - (irn) 2 ] 



E 



(18) 



This transcendental equation may be easily solved and the voltage-current characteristic (IV C) of FFO 
(due to Josephson relation the voltage is proportional to flj) flj — rj/a + fl 2 may be found, see Fig. 1, 
where results of computer simulation of Eq. (pi) and Eq. (FL8|) are presented for a — 0.2; 0.5, L = 5, T = 3. 
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On the other hand, expressing bias current rj via Qj and Q2 we arrive to exactly the same expression for 
the current-voltage characteristic derived in |l6| , Il7) : 



(2 - <5 ,n)afij(rL) 2 [l - cos(rL) cos(Trn)] 



n 2 ) 2 }[(TL) 2 - (ttti) 



212 



(19) 



So, if one needs to obtain the function Clj(r]), Eq. jl9| ) is more useful. If, however, the function Slj(r) 
(voltage versus control line current) is of importance, one can use Eq. (p"8|). When necessary (e.g., when 
r//a is of the order of unity), one can recourse to the 4-th and higher order approximations and derive 
IV C with the desired precision. 

The Eq. ( |l3| ) is the equation for slow component of the phase in the sense that it is considered in 
the steady-state limit for t — > 00. Namely such equation is required to derive different steady-state 
characteristics, such as correlation functions and spectra. Eq. ([l3]) is linear with respect to ip(x,t), 
but nonlinear with respect to 77 and T. Substituting now 77 = 770 + r]F{x,t), T = Tq + Tp{x,t) (in our 
model Tp(x,t) — arjp(x,t)) and tp(x,t) — %l>o(x,t) + ip(x,t) into ( [l3|) and linearizing it with respect to 



small fluctuations (f{x, t, 77, T) = f(x, t, 770, Fo) + 
taking into account that rd 



df(x,t,r,,r ) 
drj 



an. f (r),r ) 

drj 



V=Vo 



and r% L 



77=770 
an.j(r) ,r) 
9r 



T] F (x,t)- 



df(x,t, m> X) 



r=r, 



r F (x, *)+...), 



r=r„ 



we will get the following equation 



for the correction of the phase due to effect of fluctuations tp(x, t): 



iptt + aipt - ip xx = (r d + crr% L ) 



df(x,t) 



T) F (x,t), 



(20) 



where rd and r^ L are dimensionless dynamical resistances of FFO and control line, respectively. Following 
the linearization with respect to small fluctuations can be done if in the area of evolution of the 
fluctuating parameter 77 = tjq + r/F(x,t) bifurcation points are not located, that corresponds to the 
previously assumed condition of the smallness of fluctuations: fluctuations are so small that do not 
affect current- voltage characteristic and, therefore, do not change qualitative behavior of FFO. Equation 
( |20| ) is more general than its derivation. In spite it is obtained in the second order approximation, all 
quantities in rhs of ( poj ) are fundamental and improving the approximation up to higher order will only 
improve quantitative values of r d and r^ L and will add components with cos3£ljt, cos4f2j/; and so on 
into f(x,t). The dynamical resistances are originated from the assumption of small noise: fluctuations 
feel the system as linear if their variance is small at the scale of nonlinearity and current fluctuations are 
linearly converted into voltage fluctuations via transfer factor that is the derivative of the current-voltage 
characteristic at the working point and are not connected with " adiabatic approximation" : the noise 
in rhs of (|2C| ) is wideband, but the spectrum of "4i(x, t) depends both on properties of r]p(x,t) and the 
differential operator of the lhs of J2(i|). If we will consider the short junction limit of (|2(]) {ij) xx = 0, 
r d, L = 0)i neglect by parametric effects f(x,t) — and change variables: v = tpt> T — & r d,t, we will get 

the equation for fluctuational component of voltage v, presented in || (Eq. 4.34, page 106). 

_ 00 

The solution of Eq. ( p0|) can be expressed as a Fourier serie in space: tp(x,t) — ^2 A m (t) cos(k m x), 



where k m — irm/L. Substituting this anzats into (p0[), multiplying by cos(fc„a;) and integrating J, one 



can get the following equation for A m {t): 



d 2 A m dA m ,2 a t- i4.\ 
L a^— + k m A m = € m [t), 



dt 2 ' " dt 

where £ m (t) is the projection of the noise along the k m mode 



(21) 



2 - 60 . 



(r d + or c d L ) 



df(x,t) 



t)f(x, t) cos(k m x)dx. 



(22) 



Let us analyze correlation and spectral properties of random process £ m (t). This is nonstationary 
process due to periodic time dependence of f(x,t) and its autocorrelation function (£ m (t)£ m (t + r)) 
depends on current time t. Moreover, the product Vf{x, t) formally is a process linearly growing 

in time and does not belong to neither second nor third kind random processes that complicates the 
analysis. Here we can use one effective trick: let us analyze statistical properties of the process Cm(t) — 
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2 S ^' m f (r d + ar^ L )f(x 1 t)iip(x, t) cos{k m x)dx that belongs to the third kind if the stationary process 
tif{x, t) is white noise and later we will take derivative over Qj from the square root of the intensity of 
the noise £ m {t) . Using this procedure we can apply the standard technique |jj and obtain the correlation 
function of the second kind of the process ( m (t) : 

T* 

$c(r)= T lim o ^ J (C m {t)Cm(t + T))dt, (23) 

— T* 

that will lead to stationary delta-correlated process with some intensity, and finally we get the following 
correlation function for the process £ m ( r ) : 

$Cm(r) = a 2 {r d + or c d L f D m S{r), (24) 

where 

An = 2 ^^T ( 2 " *o,*0 {1 + H m } , (25) 



(2 - 6 , m ) 
8a 2 



dClj 



\ 



Ai n COS k n X J + I B\n COS k n X 



\i=0 



(cos k m x) 2 dx 



Let us note, that 



Dn = 



RnI 2 



P - {1 + H } 



(26) 



(27) 



for Hq — completely coincides with the dimensionalized noise intensity for a short (lumped) Josephson 
junction (see The term H m comes from down conversion of the 2-nd harmonic due to multiplication 
of cos(2f2ji + Tx) and sin(2Qjt + rx) by the noise term tif(x, t). This effect, as proven in |Q, takes place 
if the power spectral density of the process t]f(x, t) is so wide that significantly different from zero at 20, j 



and higher. The term describing the effect of the first harmonic 



dsin{n.jt+Tx) 

on., 



•qp{x,t) will not give any 



(additional) contribution, since the intensity of the fluctuational process fi(t) — siri(fljt + Tx)r/p(x,t), 
(t) = D* n 5(r)/2 will be constant as function of fij and its derivative over fij will give zero. Therefore, 
we have classical effect |j| of the additional parametric broadening of the main harmonic at the frequency 
lu = Vlj due to effect of higher harmonics (at 2£!j, and so on). 

Since for practical FFOs the second harmonic is rather weak (the output signal is nearly sinusoidal), 
and since, as it will be demonstrated below, noise components £ m (£) with m/0 have rather small effect 
on fluctuational characteristics of FFO, we for simplicity of analysis will neglect the term H m with m ^ 0. 
For m = Hq takes the form: 



[i] - Vtj/r d ) 2 



(28) 



In Fig. 2 we present the plot of the ratio between amplitudes of the second and first harmonics A 2 /Ai 
(diamonds) and the excess noise Ho (crosses) as functions of flj using the approximate Eq. ( f[8| ) for 
a = 0.5, L = 5, r = 3. It is seen, that with increase of oscillation frequency £lj the second harmonic 
amplitude becomes smaller than the first one (our approximation works better in high-frequency limit), 
but slightly increases at Fiske and Eck steps. The same qualitative behavior demonstrates the excess noise 
intensity H . Unfortunately, the expression (pq), in spite of its general form, is of restricted usefulness 
and may give incorrect results at Fiske steps. This is due to the fact that H m is derived in the second 
order approximation only, but it contains the factor l/r d . As it is seen from Fig. 1, at Fiske steps 
the approximation ( |l8| ) gives significantly underestimated values of dynamical resistance r d that leads to 
overestimated values of Ho. Contrary, at the Eck step, where approximation (18) and results of computer 
simulation nearly coincide, Eq. (^8|) gives adequate description of excess noise term Ho- In the short 
junction limit L — > one can get the following approximate expression for the excess noise term: 



Hos 



2n 2 ) 



(29) 



G 



III. FORM AND WIDTH OF SPECTRAL LINE OF FFO 



In this section we will consider the influence of noise on broadening of spectral line of the FFO. Following 
the general setup of the problem Q), we will consider the output signal at the end of the junction in the 
form: v(L, t) = ilj + v(t) + Ro cos(f2,/t + J i/(t)di), where v(t) are frequency fluctuations, and J v(t)dt is 
supposed to be slow process in comparison with cos(f2 jt) (we do not consider here amplitude fluctuations, 
Ro = const, since it is known |jj that they lead only to some noisy pedestal and do not influence the 
linewidth) . 

The correlation function of frequency fluctuations may be derived from equation ( pl| ) (see, e.g. Jlg|]). 
Since the slow component of the effective noise £ m (i) in ( pl| ) is stationary, for m^Owe can write the 
equation for autocorrelation function of A m (t), K Am [r] = {A m {t)A m {t + r)} which depends only on time 
difference r: 



d 2 K A Jr} dK Am [r] 



a- 



k 2 m K Am [r}=0, (30) 



dr 2 dr 
and should be solved with the following initial conditions: 

d2R d ^ [T] |r=0 = -a(r d + ar% L ) 2 D m /2, K Am [0] = a(r d + ar d JL ) 2 D m /(2k 2 m ). 

The correlation function of ip(L,t), if^[r], may be expressed as follows: 

oo 

K^r] = Y, K aJt}, (31) 

m— 

since cos 2 (fc m L) = 1. By the property of the correlation function, the correlation function of frequency 
fluctuations v{t) = d ^^' f ^ , ifj,[r], is the negative second derivative of if^,[r]: 

d 2 KM , , 

K M = — ^r 1 - (32) 

One can see, that for m = the correlation function K Am [r\ diverges, that reflects the only fact, that 
for the mode m = the process A$ (t) is nonstationary. Namely divergence of [r] (in our case due 
to divergence of ifA [ r ]) leads to finite linewidths of oscillators [Q: in the opposite case when K^,[t] is 
finite the linewidth will be zero. The divergence of K Ao [r] does not lead to any mathematical difficulties 

since we need to obtain (Ao(t)Ao(t + r)^ that is finite and may be derived from K Am [r], using ([32]) and 

limiting transition for m —> 0. 

Finally, one can get the following expression for the correlation function of frequency fluctuations: 



K v [t) =a(r d + ar d 1L ) 2 D x 



L p - aT j_ p—otr/2 



OO r 



-' "' - < E cos(/(q, fc m )r) - 2 f{a.k m ) sin(/(a,fc m )r) 

rn — 1 L 



(33) 



where f(a,k m ) = f(k m ) = \jk 2 n — (^) ■ If k 2 n < (^) , then one has to change sin and cos to the 
corresponding hyperbolic functions. 



For stationary and Gaussian frequency fluctuations (since Eq. (21) is linear with Gaussian noise, the 



probability distribution of A m is also Gaussian), the form of spectral line may be written as |L|] 

W v (lu) = J J exp [-x(r)] cos cordr, (34) 

— oo 

where x( T ) 1S statistical structural function that is nonnegative and even function of t: 

+r 

X(r) = \ f(T-\t\)Kv[tn- (35) 
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Substituting K v [t] ( |33|) into (|35|), we get: 

X {r) = ^r d + ar c d L fD x 

£ f [cos(/(fc m )r) + ° si g{^ )T) 



, r > 0. 



(36) 



One can check, that the statistical structural function (36) that we will use for calculation of the linewidth 
and the form of spectral line is a smooth finite function and the sum in ( |36"| ) is converging due to the 
term 

For the case of stationary Gaussian fluctuations of frequency the linewidth is defined in the following 
way H: 



An 



J exp[-x(r)]dr 



(37) 



Substituting the structural function (36) into fl37|), we get the following final expression for the linewidth: 

7T 



Ail 



J exp[-Fi(r) ~F 2 (T)]dr 
o 



(38) 



where 



F 1 (r) = -(r d + ar^) 2 D 



(39) 



F 2 (t) = l(r d + 



Ch\1 n 

ar d ) D 



2a 



oo 

-ar/2 \ 

m=l " 



cos(/(fc m )r) 



asin(/(fc m )r) 
2/(* 



(40) 



In the following we will not recourse to two known limiting cases of very fast and very slow fre- 
quency fluctuations (Lorentzian and Gaussian form of the spectral line) but perform the exact analysis 
of linewidth and spectral form on the basis of expressions (|34])-(|4"o|). 

Let us analyze the linewidth given by expressions (J3^)-(|4C|)~First, consider the case of a short Josephson 



0, f 



CL _ 



0, neglecting by parametric effects H m — 0. Then the function F 2 (t) disappears: 



junction L 

F 2 (t) = 0. If damping coefficient a is large (but, of course, rj/a > I), one can neglect the term 
(e~ aT — l)/a in (|3"§|) , (|39|) and then well-known expression for the linewidth in Lorentzian approximation 
may be obtained: 



(41) 



that in dimensional units looks like: 



A/ s — (- 
J 2 U 



R 



2 k B T 
Rn 



(42) 



It should be noted, that this expression is larger than one in Q by a factor of 7r/2 (see formula (12) in JtJ, 
we neglected by I p and Rn = Vo / I qp in their notations) . This may be explained by different definitions 
of the linewidth: we define it as the width of rectangle with the equal square Q, that for Lorentzian form 
of spectral line should give just by a factor of 7r/2 larger value than the definition of the linewidth at 1/2 
level, that was used in 0. In the following we will present all plots of the linewidth (|8|)-(fl0|) multiplied 
by the factor 2/7T for correct comparison with experiment where the linewidth is defined at the level 1/2. 

If, however, damping coefficient a is rather small, there may be significant deviation of the linewidth 
of a short junction from (41) depending on values of Do and r d (at given Do and r d the linewidth will be 
smaller), and the use of exact formula ( |38| ) for F 2 — is necessary. 

Now let us consider plots of the exact expression of FFO linewidth (|3^)-(f40|) neglecting by magnetic 
field fluctuations r^ L = 0. Our aim here is to understand how spatial modes with m/0 influences the 
linewidth, i.e. how long junction differs from the short one. The plots of the linewidth versus dynamical 
resistance are presented in Fig. 3 in dimensional units (MHz vs Ohm) for practical FFO parameters 
H: L = 76.24, a = 0.0074 (ah < 1), a = 0.04 (aL > 1), T = 4.2K, R N = 0.04 Ohm, the dynamical 
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resistance is: Rd = 0.01 — 10 Ohm; we neglected by the excess noise term Hq = 0. From Fig. 3 one can 
see, that below certain threshold curves for aL < 1 and aL > 1 have the same behavior and coincide 
with the short junction case (|42|). Above the threshold that depends on noise intensity ( p7j ) one can see 
the effect of spatial modes: the curves split and the linewidth for aL > 1 is greater than for aL < 1. 
With increase of noise intensity the threshold region located for T = 4.2K between Rd = 0.1 — 1 Ohm will 
move to the left. However, it should be noted that practical range of dynamical resistance lies from 0.001 
to 0.1 Ohm, where the effect of spatial modes can be neglected that confirms our previous assumption 
to neglect the excess noise term H m for ra/0. 

Therefore, it seems that in the practical range of parameters the linewidth of FFO may be well described 
by the following approximate expression that may be derived from (|38|)-(|40|) neglecting by spatial modes 
with m/0 and by the term [e~ aT — l)/ct: 

Mffo = \{ 2 £) + ^) 2 ^ (1 + #°)> ( 43 ) 



where Hq is given by (|2q). Now let us compare the expression (43) (multiplied by 2/n) with the exper- 



imental results [|6|], see Fig. 4. We again take for simplicity Hq = 0, as before L — 76.24, T = 4.2K, 
R N = 0.04 Ohm; for aL < 1: a = 0.0074, for aL > 1: a = 0.04. In g R% L has not been explicitly 
measured and we have used crR^ L as fitting parameter: putting the noise conversion factor a = 1 we 
have chosen R^ L to fit experimental results only at one point for Rd — > 0: for aL < 1 R^ L = 0.006 Ohm 
and for aL > 1 R^ L — 0.04 Ohm. One can see good agreement between expression ( fHf ) and the results 



of experiment. More detailed comparison of expressions ( 38 )-(^0j) and ( [43] ) with experimental results 
obtained in different layouts, taking into account the excess noise term Hq will be given elsewhere. 

In fact, Fig. 4 and the expression ([l3|) give an idea how the noise conversion factor a may be measured. 
If noise intensity is known, setting Rd <C Rd L one can get the value (aR^ L ) 2 from experimentally 
measured plot of the linewidth. On the other hand, the value R^ L is independently accessible from 
experiment. 

Let us consider the form of spectral line of the FFO. In general the spectral line consists of 
narrow and high spectral peak that finite width is originated by frequency fluctuations (nonstationary 
phase fluctuations) and broad and low pedestal due to amplitude fluctuations. If amplitude and frequency 
fluctuations are correlated but small, there will be also small asymmetric contributions both into the peak 
and the pedestal. Since in the frame of the present paper we do not consider amplitude fluctuations, 
below we will consider the form of spectral peak, that may be derived from (|34|), substituting statistical 
structural function (|36|). 

Since in practically interesting range of parameters the linewidth is well described by formula (^3|) , we 
will also neglect by spatial modes m ^ in ( fj6| ) and by (eT aT — l)/a and will get well-known expression 
for the Lorentzian form of spectral line: 

*;(„, . M 2 <^7f> , ,44) 
4tt (Af FFO /iry + lu 2 



the width of this curve at the level 1/2 is given by 2(AfpFo/' K ), where Afppo is given by (|43|). In 
Fig. 5 plots of spectral form, given by (U),(|3^) (solid lines) and the approximation ([S]) (dashed lines) 
are presented. It is seen, that for a small dynamical resistance (Rd = 0.001 Ohm, that corresponds to 
plato in Fig. 4) the curves absolutely coincide. With increase of Rd the exact expression slightly deviates 
from the Lorentzian approximation, but in all practical range of parameters (up to Rd = 0.1 Ohm), 
formula ( f44| ) gives adequate description of the form of spectral line. For larger Rd and noise intensity, 
the deviation from Lorentzian form will increase and it is necessary to use expressions (p4|) , (|3"6|) . Note, 
that in |L2| Lorentzian form of spectral line of FFO was predicted using another approach. Recently, 
Lorentzian form of spectral line of FFO has been experimentally observed in wide range of parameters 
both at Fiske and Eck steps jl3| . 

If we would separately consider the case of technical fluctuations (since these fluctuations are slow and 
narrowband the consideration may be performed in the adiabatic approximation and in this case the 
effect of additional parametric broadening of the linewidth would not appear due to the fact that the 
spectrum of technical fluctuations is much more narrow than the basic frequency of the FFO), the form 
of the spectral line will be Gaussian (for Gaussian distributed technical fluctuations) and the linewidth 
will be given as \/2tt (v 2 ), where (z/ 2 ) is variance of frequency fluctuations. It should be noted, that if 
one wishes to consider the joint effect of natural and technical fluctuations (as they coexist in real life), 
one can not recourse to the Lorentzian and Gaussian limiting cases, but the calculation of the spectral 
form should be performed using formula (p4|). 
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IV. CONCLUSIONS 



The aim of the present paper is to investigate the influence of wideband fluctuations of bias current 
and magnetic field on dynamics of FFO. We have derived analytical expressions both for the form of 
spectral line of FFO and its width. For practical range of parameters simple approximate expression 
of the linewidth, that well fits the experimental results, has been given. The appearance of excess 
noise (in comparison with the short junction linewidth) has been explained by presence of magnetic 
field fluctuations and by the so-called parametric broadening of spectral line due to influence of higher 
harmonics. It has been demonstrated that in the practical range of parameters, in the case when thermal 
fluctuations dominate, the Lorentzian form of spectral line is realized, while for larger values of dynamical 
resistance and temperature deviations from Lorentzian form may be observed. 
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FIG. 1. Current- voltage characteristic. Numerical solution of the sine-Gordon equation is presented by crosses 
and the second order approximation is given by solid line for L — 5, F = 3, a = 0.2; 0.5. 
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FIG. 3. Illustration of the effect of splitting of the linewidth for ah < 1 and aL > 1. 
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FIG. 4. Comparison of experimental and theoretical linewidths for the parameters: L = 76.24, T — 4.2K, 
Rn = O.OAOhm, (7 = 1; for ah < 1: a = 0.0074, 7?J L = 0.006 O/im; for «L > 1: a = 0.04, 7?J L = 0.04 Oftm; 
solid lines - theory, crosses and squares - experimental results. 
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